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INTEGRABLE THIRD-ORDER SYMMETRIES 
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Abstract. The complete lists of vector hyperbolic equations on the sphere and on the cone 
that have integrable third order x- and y-symmetries are presented. Integrals and Backlund 
transformations for some of these equations are found. We also find all integrable vector 
evolution equations of the third order on the cone. 



1. Introduction 

The symmetry approach to classification of integrable PDEs (see, for instance, [1-4] ) is based 
on the existence of higher infinitesimal symmetries. It is very efficient for evolution PDEs with 
one spatial variable. In particular, all integrable equations of the form 

(1-1) Ut = u xxx + $(u,u x ,u xx ), 

were listed in [5,6] (see [7] for details and proofs). For hyperbolic equations 

(1.2) U X y = ^(U,U X ,Uy) 

the symmetry approach assumes the existence of at least one x-symmetry of the form 

(1.3) u t = A(u,u x ,u xx , ...,), 
and at least one ^-symmetry 

(1.4) U T = B(u,Uy,Uyy, . . . ,). 

For example, the famous integrable sin-Gordon equation 



admits the following symmetries 



u xy = S1I1U 



1 3 , 1.8 



These equations are integrable themselves and can be found in the list of integrable equations 
of form (1.1). 
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For the general classification problem, all the functions *f?,A,B in (1.2)-(1.4) are unknown. 
They have to be found from compatibility conditions for (1.2), (1.3) and (1.2), (1.4). Such com- 
plete classification of integrable equations (1.2) turns out to be an extremely difficult problem. 
This problem is still unsolved although some partial results were obtained in [8-10]. However 
if we fix somehow functions A and B, then it is not difficult to verify whether or not the cor- 
responding function ^ exists. In particular, to find all integrable equations of the sin-Gordon 
type, one can assume [11] that both symmetries (1.3) and (1.4) are integrable equations of the 
form (1.1). 

Our goal is to apply this idea to find interesting integrable vector hyperbolic equations of 
the form 

(1.5) u xy = h u + h\U x + h 2 u y 

on the sphere and on the cone. Here u is an iV-dimensional (or even infinite-dimensional) 
vector and hi are some special scalar-valued functions. Namely these functions depend on two 
different scalar products (• , •) and (• , •) between vectors u, u x and u y . The sphere and the cone 
are defined by u 2 = 1 and u 2 = 0, respectively. Here and in the sequel for any vector a we 
denote a 2 = (a, a) and \a\ = y/a?. 

Consider equations on the sphere. The constraint u 2 = 1 imples (u, u x ) = (u, u y ) = and 
(u,u xy ) = —(u x ,u y ). From these identities it follows that equation (1.5) has the form 

(1.6) u xy = hiu x + h 2 u y - (u x , u y )u. 

Equation (1.6) is called isotropic if the coefficients involves the scalar product (• , •) only. In 
this case the hi are functions in three scalar variables (u x , u x ), (u x , u y ), (u y , u y ). It is clear 
that isotropic models are invariant with respect to the group SO(N). 

Otherwise the equation is called anisotropic. In the anisotropic case the functions hi may 
depend on six additional scalar variables (u, u), (u, u x ), . . . , (u x , u y ). All scalar products will 
be considered as independent variables (cf. [12]). As a result all our models are integrable for 
any dimension N. 

Integrable isotropic and anisotropic vector evolution equations of the form 
(1-7) u t = u xxx + f 2 u xx + fiu x + f u 

on the sphere were studied in [12, 13]. By integrability of such equations we mean the existence 
of infinite series of commuting flows 

d l u 

u Tk = g k u k + g k _i u fe _i H \- g 1 u x + g u, u { = — , 

whose coefficients ^ depend on all possible scalar products between u, ...,u k . 

The symmetry approach for vector evolution equations was developed in [12]. In particular, 
this paper contains necessary integrability conditions for equations of the form (1.7). It is 
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remarkable that these conditions are very similar to the integrability conditions (see [1,19]) for 
scalar equations (1.1). 

In [12] we have presented a complete list of isotropic equations of form (1.7) on the sphere. 
Moreover, some examples of integrable anisotropic equations were found there. One of these 
examples is the following equation [14]: 

(1.8) u t = (u xx + -(u x , u x )uj +^(u,u)u x , u 2 = l. 

Here (a,b) = (a, Rb), where R is an arbitrary constant symmetric matrix. One can assume 
that R = diag(ri, . . . , r^). Equation (1.8) has a Lax pair whose spectral parameter lies on an 
algebraic curve of genus 1 + (N — 3)2 N ~ 2 . If N — 3, then (1.8) is a commuting flow for the 
famous Landau-Lifshitz equation. 

A complete list of anisotropic integrable equations of the form (1.7) on the sphere was pre- 
sented in [13]. Equations without any constraints on the length of the vector u were investigated 
in [15-17]. Neither evolutionary nor hyperbolic integrable equations on the cone have never 
been considered. 

In Section 2 we find all hyperbolic equations (1.5) that have symmetries from the lists of the 
papers [12, 13]. 

Section 3 is devoted to integrable vector equations on the cone u 2 = 0. In Proposition 2 we 
describe all non-triangular integrable vector evolution equations of the form (1.7) on the cone. 
We also find hyperbolic equations (1.5) that have both x and y integrable symmetries of the 
third order from the list of Proposition 2. 

Acknowledgments. The authors are grateful to V. Adler, and R. Yamilov for useful discus- 
sions. The research was partially supported by the RFBR grant 11-01-00341-a. 

2. Integrable hyperbolic equations on the sphere 

Two equations (1.5) on the sphere u 2 = 1 are said to be equivalent if they are related by a 
composition of the following point transformations: 

1. The scaling transformations of the form: 

(2.1) x' = ax, y' = (3y; 

2. The interchanging y -H- x; 

3. The transformations of the second scalar product 

(2.2) (u,v) ->■ a{u,v) + p(u,v), 

where a and j3 are arbitrary constants. In general the constants in the equivalence transforma- 
tions may be complex. 
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Proposition 1. If equation (1.6) has integrable symmetries of the form 

(2.3) u, = u xxx + f 2 u xx + fiu x + f u 

and 

(2.4) U T = Uyyy + Q^Uyy + Q^U y + Q Q U 

then it is equivalent to one of the following equations: 

(2.5) u xy = yjl - u 2 yS Jl + \u x \- 2 u x - (u x , u y )u, 

(2.6) = ((u,Uy) + \Jl + (u,u)\u x \- 2 ip) - (u x ,u y )u, 

Lp = yj (U, U y ) 2 + (U, U){\ ~ (Uy, Uy)), 

(2.7) u xy = \j\-u\ Uy - (u x , Uy)U, 

(2.8) u xy = Uy ^ u x + J (u, u) - u 2 x Uy - (u x , u y )u, 

[u, u) 



(u, u x ) + J (u, u x ) 2 + («, u) (l - (u x , u x )) 

(2.9) U X y = r Uy - (U X , Uy)U, 

\U, U) 

(U, Uy) (U, U X ) + if) 

(2.10) U X y = ~U X H r Uy - (U X , Uy)U, 

(it, It) \U,U) 



lj>= yj{u, U x ) 2 + (it, U) 2 - (u, U)(U X , U x ) . 

For the sake of brevity we do not present here the symmetries for equations (2.5)-(2.10). 
These symmetries can be easily calculated with the help of special computer software (see 
for example [20]). Another way is to try one by one the evolution equations from [12, 13] as 
symmetries. 

Remark 1. Any anisotropic equation admits the reduction (•,•)—>■ 7( • , •) (isotropic limit). 
Equation (2.6) is reduced to (2.5) and the remaining anisotropic equations shrink into (2.7). 

It is interesting to note that symmetries for some anisotropic equations are isotropic. In 
particular, equation (2.9) has the following isotropic y-symmetry 



2 



{Uy, Uyy) 3 U yy 



o v-«; -yy, i gg, ? . 

yyy ,.2 w ' 9 ,.2 y 



(the vector Schwartz-KdV equation [18]). 

2.1. One-component reductions. All equations (2.5)-(2.10) are vector generalizations of 
different special cases for known scalar integrable equation [10,21] 



(2.11) u xy = P(u)yJu x + l y Jul + a. 
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Here P is any solution of differential equation P' 2 = AiP 4 + A2-P 2 + A3, where Aj are arbitrary 
constants. The generic solution of this PDE is the elliptic Jacobi sine. 

Indeed, consider the case N = 2. Let the scalar products be u 2 = u\ + u 2 , and (u, u) = 
au\ + bul, a 7^ b. Using the trigonometric parametrization u = {cos w, sin w} of S 1 , one can 
rewrite equations (2.5)-(2.10) in terms of the scalar variable w. 

This scalar reduction for (2.5) is the equation w xy = yjl — w 2 + w 2 equivalent to (2.11) 
with P — 1. The reduction of equation (2.7) is given by a/1 - w*. 

The scalar reductions for the anisotropic equations (2.6), (2.8) and (2.9) can be brought to 
the form (2.11) by the following point transformation 

dw f df 



(2.12) . , , 

J V® cos2w + f3 J ^/(l + / 2 )(a + 6/ 2 ) 

where a = a + (3, b = (3 — a and / = tan w. 

The reductions for equations (2.6), (2.8) and (2.9) are given by 

(2.13) v xy = : J^^l + v 2 x ^l-abv 2 , 



(2.14) v xy = v y w'{v)y/l-vl, 

(2.15) v xy = ^^l-abv 2 , 

respectively. Equations (2.14) and (2.15) can be reduced to (2.11), where a = 0. For equation 
(2.10) we have 



(2.16) v xy — Vy^/Y— abv 2 , where w'(v) = a cos2w + f3. 

Note that (2.13) is the generic equation of the form (2.11) and therefore (2.6) seems to be 
the most interesting integrable model found in Proposition 1. 

2.2. Two types of integrability. It is known that there exist two different classes of inte- 
grable equations (1.2): the Liouville and the sin-Gordon type equations. Equations of the first 
class (also named as Darboux integrable equations) have so-called x and ^/-integrals. A function 
/ depending on u, u y , u yy , . . . , u y n ^ is said to be an n-th order x-integral of hyperbolic equation 
(1.2) if the identity 

dx 

is valid for any solution u(x,y) of (1.2). For example, the function / 

integral for the celebrated Liouville equation u xy = exp(it). The sin-Gordon type equations are 
integrable by the inverse scattering method (see for example [22]). They have no local x- or 
^/-integrals. 

For vector equations of the form (1.5) the x-integrals are functions I depending on several 
scalar products of vectors u, u y , u yy , .... Equations (2.7)-(2.10) possess integrals. For instance, 
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the function 



Jl = (u x - Uy/l - U 2 X , U x - Uy/l-U^, 

is the y-integral for equation (2.7); equation (2.10) has the following y-integral: 



J 2 = #±^, where h = ^A±± , 




and ip is presented in (2.10), and so on. It is remarkable that isotropic equation (2.7) has the 
anisotropic integral. 

Equations (2.5) and (2.6) have no integrals. 

2.3. Backlund transformations. The first order auto-Backlund transformation (BT) for the 
hyperbolic equation (1.5) is a pair of the ordinary differential equations of the form 

(2.17) u x = fiu + f 2 v + f 3 v x , u y = g x u + g 2 v + g 3 v y , 

compatible with equations (1.5) and 



with the same functions hi as in (1.5), where one must change u in v. The functions /$ in 
(2.17) depends on various scalar products of the vectors u, v and v x , whereas gi are functions 
in various scalar products of the vectors u, v and v y . 

If a solution u(x,y) of equation (1.5) is given, then one can find another solution v(x,y) of 
the same equation by solving ODEs (2.17). Suppose that the coefficients of (2.17) depend on 
arbitrary "spectral" parameter A. Then starting from a simple solution of (1.5) and applying 
(2.17) several times, we can construct a multi-parametric family of solutions. The so-called 
multi-solitonic solutions for classical integrable systems can be found in this way. By that 
reason the existence of BT with spectral parameter may be regarded as an indication of the 
complete integrability. 

To justify the integrability of equations (2.5) and (2.6) we present for them Backlund transfor- 
mations with the spectral parameter. Equation (2.5) has the following BT: 



v xy = h v + hiv x + h 2 v. 



y 



u 



'x 



f [(«, v x )(u + v)- ((«, v) + l)v x ] , 



(2.18) 



u 




((«,«) + 1) (A + 1) 



X(U,Vy) 



- g \v - v 



where 




1 1 1 l ((u,v) + l + 2X)(v 2 x + l) 

A + {u,v) + l + \\l vl((u,v) + l) 
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The Backlund transformation for equation (2.6) has the following form 

q + fh [ (u,v x ) 

U x =-, r- 7 —-{U + V)-V x , 

(2.19) i(u,v) + l J 

u y = («, v y )u + g((u, v)u - v) - Vy, 

where 

q = A((w, u) + 1) - (u,v), f = y/l + (v,v)K|- 2 , h = ^/q^^u/u^v^ , 

3. INTEGRABLE EQUATIONS ON THE CONE 

V. Adler has put our attention to the fact that integrable hyperbolic equations on the cone 
u 2 = cannot be drawn from integrable equations on the sphere u 2 = r 2 by the limit r — > 0. 
Indeed, applying the transformation u — > r~ x u to equation (1.6), we see that this limit does 
not exist because of the term (u x , u y )u in (1.6). In this section we consider integrable evolution 
and hyperbolic vector equations on the cone. 

The constraint u 2 = implies (u, u x ) = 0, (u,u xx ) = —(u x ,u x ), (u,u xxx ) = —3(u x ,u xx ) 
and so on. It follows from these relations that any third order vector evolution equation on the 
cone has the following form 

(3.1) u t = u xxx - 3 Uxj ™ xx u xx + /itt a + f u, 

where /o and f\ are scalar- valued functions depending on three isotropic variables u 2 , (u x , u xx ), 
u 2 x and six anisotropic variables (u, u), (u,u x ), . . . , (u xx , u xx ). 

The class of equation (3.1) admits the following point transformations. 

1. The scaling transformations: 

(3.2) x' = ax, t' = a 3 t; 

2. The Galilean transformation: 

(3.3) x— > x + ct, u t u t — cu x ; 

3. The linear transformations of the scalar products: 

(3.4) (u, v) — > a(u, v) + b(u, v), (it, v) — >■ c (u, v), 

where a or b may be zero, c ^ 0; 

4. The transformations of the form: 



(3.5) 

where / is an arbitrary function. 



u = vf((v,v)), 
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5. If equation (3.1) is homogeneous: f k (Xu, Xu x , \u xx ) = f k (u,u x ,u xx ), VA, k — 0, 1, 
then it admits the transformation 

(3.6) u = v exp(a x + bt). 

If two equations of the form (3.1) are related by some of transformations (3.2)-(3.6) then these 
equations are said to be equivalent. 

There exist so called triangular equations on the cone. 
Example 1. The following equation 



... -\ n { u x, u xx) . 3 It 

(3.7) u t = u xxx - 3 - u xx + -— + fou 

ui 2 ui 



on the cone, where the fo is arbitrary, can be decomposed into the same equation with fo = 
for a vector variable w and the additional equation 

f(w x , w xx ) n r x \ 3 fw 2 (w x ,w xx )r x r 2 .\ 
r t = r xxx - 3r xx 1 2 xxJ +2-l)+-r x (^f + 4 l 2 xx) - + 4-§ + r / 

for a scalar variable r. Here it = r (io + c) , where w 2 = 1 and c is a constant vector such that 
c 2 = —1 and (w,c) = 0. The equation for w is integrable on the sphere [12] but the whole 
system is not integrable for the generic function f in the equation for r. Equation (3.7) has 
infinitely many conservation laws that can be expessed in terms of w and its derivatives. While 
the equation with the generic function f Q has no higher symmetries. 

Proposition 2. If a non-triangular equation (3.1) possesses an infinite series of symmetries, 
then it is equivalent to one of the following equations: 

, ^ Au x ,u xx ) 3 (u 2 (u x ,u xx ) 2 \ , / 7T 

3.8 u, = u xxx - 3 l x \ xx, u xx + - + \ x \ xx) u x + by/l + aulu, 

v 2 2\u 2 x u x (l + au 2 x )J 



(u x ,u xx _ / u 



,2 



(3.9) u, = u xxx - 3 V - xx u xx + 3 -^f + 3bu x ) u x - 3bu 2 x (u x , u xx )u, 

(3.10) u t = u xxx - S^^-u xx + I ( ™f + K '^ )2 + 

< u% u 2 x 



/Q -I -I \ Q U XX) 

\o.ll) U t — U xxx — 6 5 U xx -\- 



| 3 (u 2 xx | (n^-u^) 2 _ ((n a ,Mxg) + (vU) K,-»x) ^ 
where o, 6 and c are parameters. 

Remark 2. Equations (3.8) with 6 = 0, (3.9) with 6 = 0, (3.10) and (3.11) are also integrable 
on the sphere. 
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Proposition 3. If a hyperbolic equation of the form (1.5) on the cone has x- and y-symmetries 
from the list of Proposition 2, then it is equivalent to one of the following equations: 

(3.12) u xy = u 2 x u y , 

(3.13) u xy = u 2 x sja+ \u y \- 2 Uy, 

(3.14) u xy = u 2 y u x + u x u y + u 2 x u 2 y u, 
where a is a constant. 

All equations (3.12)-(3.14) have integrals. For example, equation (3.14) has the following 
vector integrals: 

d u x -u 2 x u _ d u y -u 2 y u _ q 

dy \u x \ dx \u y \ 

Integrating these equations we reduce (3.14) to a pair of the ordinary equations: 

(3.15) u x = u 2 x u + \u x \v(x), u y = u 2 u + \u y \w(y), 

where v and w are arbitrary vector- functions such that v 2 — 1, w 2 — 1. 

Equations (3.12) and (3.13) have the following y-integral / = u\. This allows us to find the 
general solutions of the equations in quadratures. 
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